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PROJECT  BRISHAM  wis  a  Department  of  Defease  data 
collection  effort  wherein  890  kHz  wile  samples  of  the 
HF  spectrum  were  aogiired  using  a  25  Hz  sampling  rate 
with  2.8  kHz  resolution  for  a  period  of  up  to  2.4 
minutes.  The  main  tacust  of  this  paper  has  been  the 
visual  examination  of  the  data  on  a  graphics  display 
with  the  idea  of  identifying  drastic  and  unexpected 
changes  in  either  individual  signals  or  the  entire 
wavefront;  those  changes  being  iue,  at  least  in  part, 
to  signal  transform tion  resulting  from  various 
ionospheric  phenomena.  In  addition,  an  in-depth 
statistical  study  was  conducted  on  several  of  the 
specific;  signals  witn  the  hope  of  aiding  in  the 
formulation  of  an  algorithm  to  efficiently  automate 
the  receiving  system1 s  sensitivity  settings. 
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I.   OBJECTIVES 


The  objectives  of  the  project  were;  (1)  to  continue  the 
visual,  or  graphical,  analysis  of  the  BRI3fl&a  data  (the 
initial  effort  of  which  -fas  reported  in  Ref.  1)  in  an 
attempt  to  verify  anomalous  signal  behavior  already 
documented  and  to  discovec  any  new  anomalies  whica  might  be 
present;  (2)  toward  that  end,  to  make  whatever  changes  were 
necessary  in  both  the  programs  and  data  nandling  procedures 
to  meet  either  new  systen  constraints  or  to  enhance  over-all 
project  efficiency;  and  (3)  to  analyze  the  data  from  a 
statistical  standpoint  (via  two  entirely  new  computer 
programs)  in  an  effort  to  determine  how  often  and  by  what 
amount  (in  decibels)  several  different  classes  of  signals 
change. 


The  format  of  the  lata,  particularly  that  of  the  2. '4 
minute  time  duration  limitation  (although  there  ware  several 
time-wise  contiguous  data  sets),  aade  long  tern  analysis, 
both  visual  and  statistical,  impossible.  In  addition,  the 
daily  intercept  times  for  the  data  sets  precluded  any 
significant  effort  beiag  made  to  correlate  signal 
characteristics  from  even  approximately  the  same  time  of  the 
day. 


II. 


INTRODUCTION 


It  is  well  documental  that  as  radio  signals  traverse  the 
ionosphere  (including  signals  bouncad  off  of  one  of  its 
layers)  they  ara  often  cnanged  in  a  Banner  that  goes  beyond 
simple  attenuation  due  to  the  distance  travelled.  Both  the 
make-up  of  the  ionosphere  itself  (hotf  many  ions  ire  present 
at  any  one  tine,  the  elusive  5poradic-E  phenomenon, 
night-day  asymmetry  and  the  geomagnetic  non-reciprocity  to 
name  but  a  few)  and  the  more,  dramatic  disraptions  of 
ionospheric  behavior  (solar  flare  caused  magnetic  storm 
activity,  polar  cap  absorption  and  tae  auroral  displays  are 
three)  team  with  the  aore  common,  every  day  paenomena  of 
"Faraday  Rotation"  and  noise-producing  lightning  storms  co 
alter  any  signal  with  which  any  or  all  of  taem  come  in 
contact.  The  results  can  De  severe  attenuation,  pronounced 
amplification,  "ducting",  "whistlers",  and  fading  (sometiaes 
called  the  Dellinger  3ff ect) .  References  2-5  were 
excellent  sources  for  ionospheric  phenomena  inforaation. 

It  was  felt  that  the  application  of  graphical  techniques 
to  signals  which  had  traversed  the  ionosphere  night  add  some 
insight  into  the  problem  of  exactly  what  happens  to  these 
signals. 


For  the  most  part,  the  visual  presentation  techniques, 
shown  in  various  stages  in  Fig.  1,  used  in  this  study  ware 
the  same  as  those  presented  in  Ref.  1.  However,  some 
changes  have  been  made  and  will  be  discussed  in  detail  in 
the  appropriate  section.  ks  in  Ref.  1,  interesting  sets  of 
scans  have  been  reproduced  as  graph  plots  (using  the 
VEBSAIEC  plotter)  for  further  analysis  and  presentation  in 
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Figure   1  -   ANOMALY  DATA  PROCESSING  STAGES 
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this  paper.   All  thirty-four  of   the   previously   unanalyzed 
complete   data   sets   wars   studied   with  varyiag  degrees  of 
success  concerning  anomalies  (or,  least,  perceivel  to  be  so> 
found  therin. 

The  swept-tuned  receiver  and  associated  equipment  used 
to  acquire  and  record  the  data  sets  is  an  operational 
receiver  equipment  used  foe  other  than  experimental  data 
gathering.  Its  proper  operation  requires  frequent  tuning  of 
numerous  sensitivity  adjustments.  Ho*  often  and  by  what 
amount  these  attenuators  and  thresholds  are  changed  have  a 
profound  effect  on  the  performance  of  the  equipment. 
Therefore,  programs  MIN/SAX  1  and  MIN/M&X  2  were  written  in 
an  effort  to  accumulate  eaough  statistical  data  from  the 
sets  to  warrant  the  recommendation  of  some  appropriate 
algorithm  for  controlling  the  attenuator  adjustments. 
Section   IV  of  this  paper  will  address  itself  to  that  topic. 


A.   NATORE  OF  THE  DATA 


The  data  was  gathered  using  a  swapt-tuned  receiver  with 
either  an  omnidirectional  or  12  degree  beam  antenna.  The 
data  sets  were  produced  oy  sampling  the  receiver  output  at  a 
25  Hz  sampling  rate  as  it  scanned  downward  through  an  390 
kHz  band.  Each  data  set  consists  of,  at  most,  3597  scans 
(each  scan  being  composed  of  318  contiguous  2.3  kHz  wide 
bins  moving  down  in  frequency  as  the  bin  number  increases) 
placed  on  7194  cards.  Cird  one  gives  amplitude  ii formation, 
in  2  dB  high  quantization  levels,  for  bins  1  througn  150 
while  information  on  bins  161  through  318  is  found  on  card 
two.  The  quantization  levels  are  given  in  5-oit  binary 
words  with  even  numbered  bin  information  found  in  rows  11 
through  3  while  odd  numbered  bin  information  is  in  rows  5 
through   9.    In  addition,  each  card  has  the  card  number  and 
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set  number  punched  in  row  12  foe  accounting  purposes.  With 
each  set  there  is  a'  paraneter  card  which  delineates  the  date 
and  time  of  signal  acquisition,  base  frequency,  ajitenna 
configuration  and  attenuation  applied. 

Of  approximately  150  data  sets  recorded,  forty-four  were 
made  available  for  this  study.  Figure  2  is  a  table  of 
pertinent  information  concerning  the  forty-two  conplete  sets 
analyzed  in  Ref,  1  and  this  paper. 


B.   CHANGES  IN  THE  GRAPHICAL  ANALYSIS  PROCEDURE 


Equipment  reconfigurations  since  the  graphical.  analysis 
procedure  was  formulated  aad  reported  in  Ref.  1  and  efforts 
to  make  the  data  .handling  more  efficient  demanded  several 
changes  both  in  procelure  aad  in  the  two  nai.Q  analysis 
programs,  Anomaly  A  and  Anomaly  6.  Specific  information 
concerning  the  workings  of  the  two  programs  may  be  found  ia 
Ref.  1,  particularly  in  section  II  B. 

First  of  all,  the  sheer  number  of  binary  cards  to  be 
transcribed  onto  magnetic  tape  (nearly  a  quarter  of  a 
million)  obviated  the  card-to-tape  techalgue  used 
previously.  A  much  faster  and  far  nre  reliaoLe  procedure 
was  initiated  wherein  tie  IBM  360/67  system,  belonging  to 
the  w.  R.  Church  Computer  Center,  was  used  after  normal 
working  hours.  The  cards  were  first  read  onto  9-track  tapes 
(800  BPI,  even  parity)  aad  then  transcribed  to  the  7-track 
tapes  (556  BPI,  odd  parity)  required  for  the  XDS-9300 
computer.  Partial  dumps  were  made  of  all  sets  from  both  the 
9-  and  7-track  tapes  to  ensure  that  the  cards  read  onto  the 
tapes  properly  and  that  tae  tape-to-tape  translation  was 
successful.  The  programs  used  to  read  the  cards  Dnto  tapes, 
translate  the  9-  to  7-tract  and  dump  the  tapes  for 
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Figure   2  -   DATA  SET  INFORMATION 
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quality  control  can  be  found  at  the  ead  of  this  paper. 

This   new   system   o£   lata   loading   resulting   in  the 

requirement  to  change  the  sections   of   the   Anomaly   A  and 

Anomaly   B   programs   referring   to   data   read-in.   This  was 
accomplished  with  little  difficulty. 

Next,  a  metasymbol  subprogram  was  added  to  sach  of  the 
analysis  programs  to  enable  entire  files  (data  sets)  to  be 
skipped  when  a  spacific  file  reguirei  was  located  other  than 
first  on  a  tape.  By  sailing  the  routine  "FORSCN"  a 
keyboard-entered  number  of  files  (NPILE  =  )  wouli  be  passed 
over.  This  proved  to  be  extremely  valuaole  as  a  time  and 
labor  saving  device.  The  routine  "BAK5CN"  was  also  available 
for  backing  over  entire  files  but  cnronic  tape  drive 
malfunctions  negated  the  use  of  this  even  aore  valuable 
tool. 


One  final  major  program  change  was  dicta 
computer  center  equipment  reconfiguration  tfhs 
CALCOMP  563  plotter  was  no  longer  available.  The 
to  be  used  was  the  PDP-11  computer  in  conjunctiD 
VERSATEC  plotter.  This  caange  required  modificati 
analysis  program  sections  dealing  with  the  process 
output  data  (desired  plots  from  the  &GI-10  screen) 
the  very  end  of  the  main  program,  and  the  generat 
new  subroutine  ,  called  the  Header  Subroutine,  w 
the  header  data  into  the  proper  format  for  the 
unit.  In  addition,  a  PDP-11  plot  package  (used  to 
the  data  £ape  with  the  plotter)  had  to  be  e 
modified  to  meet  the  demands  of  this  anoaaly  data, 
of  the  revised  Anomaiy  &  and  B  programs,  the 
program  for  calling  "FORSCN"  and  "BAKSCN"  and 
plot  routine  (called  "aa.c"  foe  Anomaly  A  plot 
routine  alone  is  listed  since  the  one  utiliz 
Anomaly   B   graph   plots   differs   only   in   a   fa 
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approximately  UDO  lines  comprising  "aa.  c")  can  De  found  at 
the  back  of  this  paper.  Sibroutines  3INP,  SINPUT  ,  FNS  and 
VCD  were  not  changed  aad,  hence,  are  not  listed  in  this 
paper  but  may  be  found  ia  Ref.  1. 


C.   ANOMALY  A  AMD  B  PLOTS 


The  plots  which  have  resulted  from  tha  graphical 
analysis  portion  of  tais  study  require  some  prior 
explanation.  Dua  to  tha  data  handling  restrictions  imposed 
by  the  XDS-9300  computer,  only  a  small  portion  of  tha  data 
in  each  set  could  be  viewed  at  any  one  time.  This  resulted 
in  some  difficulty  in  visualizing  what  could  be  long-term 
(time  wise)  or  wide  (freguency  wise)  phenomena  and  demanded 
that  each  set  be  analyzed  several  times  before  the  entire 
set  had  been  viewed. 

The  Anomaly  A  plots  consist  of  twenty  scans  (signal 
wavefronts  sampled  0.D4  seconds  apart  -  or  further  if  scans 
were  deliberately  skipped  to  view  long-term  effacts)  with 
the  most  recent  scan  ia  time  falling  at  the  bottom  of  the 
picture.  The  scans  could  be  brought  onto  the  scraen  in  any 
number  up  to,  and  including,  twenty  at  a  time.  ^s  shown  in 
Fig.  3,  the  vertical  axis  is  marked  off  in  5  d8  increments 
which  apply  only  to  tas  bottom  scan  (called  "base  scan"). 
The  horizontal  axis  is  divided  into  ten  23  kHz  wide  sections 
each  composed  of  ten  fregjency  bins  2.8  kHz  wide  for  a  total 
of  100  frequency  bins  and  230  kHz.-  This  width  can  be  varied 
to  other  than  100  bins  if  desired.  The  leftmost  Din  (termed 
the  "base  bin")  is  selected  by  keyoorad  input,  thus  allowing 
the  analyst  to  slide  up  and  down  tha  data  set 
(frequency-wise)  at  will.  The  scale  along  the  horizontal 
axis,  then,  consists  of  iacrements  of  ten  bins  starting  with 
the  "base  bin"  number  at  the  origin.   The  term   "oase   freg" 
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refers  to  the  base  freguency  of  the  data  set  and  decreases 
in  2.8  kHz  increments  as  the  bii  number  increases.  Shading 
has  been  provided  on  some  of  the  plots  to  highlight  the 
suggested  anomaly. 

The  Anomaly  B  plots  are  a  combination  of  the  Anomaly  A 
plot  technique  and  a  modification  which  enables  the  analyst 
to  view  selected  bins  (signals)  in  the  time  vs.  amplitude 
domain.  At  the  bottom  of  each  plot  are  five  scans  (Anomaly 
A  style)  used  to  aid  in  the  selection  of  the  bins  to  be 
plotted  in  the  upper  traces  (see  Fig.  4).  In  addition  to 
the  restricted  number  of  scans  to  be  seen  at  one  time  (only 
5)/  computer  Limitations  required  that  the  number  of 
frequency  bins  presented  be  reduce!  to  eignty.  Other  than 
those  changes,  discussion  of  the  Anomaly  A  plots  is  also 
valid  for  the  lower  set  of  Anomaly  3  plots.  The  upper  set 
of  Anomaly  B  plots  is  a  collection  of  five  keyboard-selected 
signals  chosen  from  among  the  scans  below.  The  traces  are 
numbered  one  to  five,  from  top  to  bottom  respectively.  The 
horizontal  axis  is  a  time  line  composed  of  eighty  data 
points  at  least  0.04  seconds  apart  (depending  upon  the 
number  of  scnas  skipped,  if  any)  with  the  oldest  point  at 
the  extreme  right  side  of  the  trace)  with  the  hasi  marks  on 
the  vertical  axis  refering  to  10  dB  amplitude  levels. 

As  se.en  to  some  extent  in  the  Anomaly  A  plots  and,   more 

markedly,   in   the  Anomaly  B  plots,  the  scans  loot  fractured 

in  appearance.   This   was   not   the   result   of   intentional 

programming,    rather    it    was   the   product   of   the   new 

'    PDP-11/VilRSATEC    plot   routine    conoina tion. 
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III. 


GRAPHICAL  ANALYSIS 


A.   DISCUSSION 


Prior  to  initiation  of  the  graphical  analysis  process 
via  programs  Anomaly  A  ana"  3,  an  enumeration  of  many  of  the 
possible  "anomalies"  expected  was  made.  This  was  not  to 
imply  either  that  signals  which  had  Deea  markedly  affected 
by  the  ionosphere  would  only  exhibit  these  graphical 
manifestations  of  anomalous  performance  or  that  all  forms  of 
signal  anomalies  would  be  seen  by  the  analyst.  Tais 
enumeration  procelure  was  merely  an  attempt  to  'visualize" 
what  effects  those  siga al-changiag  ionospheric  phenomena 
discussed  in  the  introduction  might  have  on  the  lata  sets 
yet  to  be  analyzed.  Although  no  attempt  was  made  to  make 
this  study's  "anomalies"  conform  with  those  iieatified  in 
Ref.  1,  there  has  been  aa  attempt  to  verify  any  previously 
seen. 

The  graphical  effects  on  the  signal  waveforns  expected 
were: 

1.  a  sudden  enhancenea t/depression  of  one  oc  more 
signals  (long  or  short  term) ; 

2.  an    abrupt   discontinuity   across   an   entire 
wavefront  or  portion  thereof; 

3.  a   signal   "bum?"  or  "valley"  moving  either  up 
or   down   in   frequency    as    the    wavefront 
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progresses  (Fig.  4 ,  Ref.  1  refers)  ; 

4.  a  shift  in  frequency  of  one  or  mora  signals 
(either  short  tera,  long  teem  or  aperiodic  in 
nature)  ; 

5.  an  unex? lainable  enhancement  of  one  or  more 
signals  with  a  concomitant  depression  of  one 
or  more  signals  at  the  same  time  in  the  same 
data  set; 

6.  a  complete  loss  of  continuity  of  the  signal 
wavefront  for  variable  periods  of  time; 

7.  unexpected  variations  in  one  or  mora  signal 
contours  as  the  wavefront  progresses; 

8.  any  out-of-the-ordinary  graphical  implication 
of  anomalous  performance. 

Effects  1,  2,  4,  5,  and  8  were  discerned  in  the  present 
study  while  several  of  these,  plus  effects  3,  5,  and  7 ,  were 
reported  in  Ref.  1 . 


B.   ANOMALY  A  ANOMALIES 


As  shown  in  Ref.  1  aid  verified  in  this  paper,  most  of 
the  data  sets  (when  viewed  on  the  graphics  terminal) 
exhibited  a  great  measure  of  uniformity  in  signal  makeup  in 
the  time  domain  and  did  not  yield  any  obvious  graphical 
anomalies.  In  fact,  fulLy  fifty  per  cent  of  tae  BRIGtiAM 
data  sets  failed  to  provide  any  repoctaDle  irregularities  in 
signal  pattern  with  the  other  fifty  per  cent  giving  fairly 
uniform   signal   wavefronts   for   most  of  the  time  recorded. 
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The  latter  group  did,  however,  exhibit  anomalous  (or,  at 
least,  assumed  to  be  so)  characteristics  which  were 
verifications  of  performaice  reports!  in  Ref.  1  or  were 
newly  discovered.  The  uniformity  mentioned  above  ailed 
immeasureably  in  the  detection  of  anomalies  sinoe  even  a 
slight  wavefront  variation  could  be  discerned  ii  the  midst 
of  otherwise  uniform  signal  progression. 

Figures  5  and  5  are  two  examples  of  minor 
discontinuities  (affecting  30  to  40  bins  -  approximately  100 
kHz)  which  appear  to  be  the  result  of  a  multi-signal 
amplitude  enhancement  for  a  very  short  period  of  time.  This 
effect  was  found  scattered  through  almost  half  of  the 
thirty-four  data  sets  viswed  via  Anomaly  A.  Figure  7  is  an 
example  of  aperiodic  discontinuities  which  ran  across  almost 
all  of  the  entire  data  set  bandwidth.  This  changs  did  not 
occur  over  the  entire  time  duration  of  the  set  and,  as 
mentioned  above,  could  not  be  traced  over  its  entire 
frequency  spectrum.  This  type  of  discontinuity  was  also 
found  in  data  sets  33  and  80.  Both  of  the  discontinuity 
examples  were  anticipated  although  the  latter  fora  holds  the 
most  promise,  anomaly-wise. 

Fading  was  another  effect  expected  but  tas  examples 
found  were  somewhat  surprising  in  that  it  was  anticipated 
that  the  entire  fregueicy  bandwidth  of  a  signaL  wavefront 
(890  kHz)  would  be  affected.  This  pcoved  not  to  be  the 
case,  with  all  of  the  examples  of  fading  implicating  at  most 
20  frequency  bins  at  a  tine.  Figure  3  shows  the  signal 
centered  on  bin  148  fading  in  and  out,  a  process  which 
lasted  several  seconds.  Figure  9  shows  a  murh  smaller 
time-duration  fading  centered  on  bin  24. 

A  third  typs  of  possioLe  anomalous  activity  is  that  of 
the  shifting,  up  or  down  in  frequency,  of  eithar  a  single 
signal  or  entirs  390  kHz    wide  signal  scans.   The  precise 
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Figure   8  -   FADING  EXAMPLE  §1 
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explanation  of  this  phenomenon  is  not  available  but  it  must 
be  kept  in  mini  that  naif un:tions  in  the  receiver  system 
itself  could  be  the  cause  rather  tian  some 
ionosphere-related  activity.  Figure  10  shows  a  possible 
left  shift  (approximate! p  3  bins)  and  distortion  of  two 
side-by-side  signals.  Figures  11  ana  12  are  examples  of 
apparent  frequency  shifts  of  160  bins  both  of  wnich  turned 
out  to  be  data  set  problems.  Since  each  scan  is  composed  of 
two  data  cards  of  binary  coded  amplitudes  (bins  1  -  160  on 
card  one  and  bins  161  -  318  on  card  two),  the  ordar  in  which 
the  two  cards  appear  and,  as  importantly,  the  existence  of 
both  cards  for  aach  scan  ara  absolutely  vital  to  the  correct 
representation  of  the  signal  wavefront.  For  example,  if 
card  one  of  a  scan  is  missing,  than  card  one  of  the  next 
scan  is  taken  as  card  two  of  the  formar.  This  has  a  domino 
effect  which  will  graphically  manifest  itself  as  an  apparent 
shift  of  160  bins  in  all  subsequent  scans.  Figure  11 
illustrates  precisely  tais  effect.  A  manual  search  of  the 
missing  from  the  set.  Figure  12  is  an  interesting  twist  on 
this  problem  ia  that  the  signal  apparently  shifts!  160  bins 
to  the  right  (or,  down  ia  freguancy)  and  then,  three  scans 
later,  shifted  right  bac<  again.  another  search  of  the 
cards  showed  that  ,as  expected,  two  cards  were  missing,  one 
to  cause  the  initial  160  bin  right  shift  and  a  second  to 
return  the  data  cards  to  their  correct  sequence  three  scans 
later.  One  last  exampLe  of  signal  shifting  is  by  far  the 
most  interesting  of  the  lot.  As  illustrated  by  Fig.  13, 
there  is  a  shift  down  in  freguency  of  approximately  12  bins 
(34  kHz)  affecting  the  entire  318  bins  (890  kHz  bandwidth) 
of  aper.iodicaliy  spaced  scans.  This  phenomenon  was  saan 
over  the  entire  length  of  the  data  set  at  apparently  random 
intervals.  Data  sets  138  and  143  also  exhibited  this 
aperiodic  frequency  shifting,  although  tae  number  of  bins  by 
which  the  scans  were  shifted  varied  slightly  over  the  three 
sets. 
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28 


u 

"Z 

w 

o* 
w 

Pi 


c 


*<v 


<ggpp>m< 

~  -  -Mi 


R 


H 

c 

•H 

pq 

cu 

w 

cd 

PQ 

A 

CN 

N 

:»fc 

PC 

S 

w 

hJ 

00 

eu 

H 

S 

<d 

cr 

X! 

cu 

w 

^ 

fa 

o 

z 

0) 

H 

en 

H 

CO 

fa 

pq 

H 

33 

A 

CO 

o 

H 

tH 

1 

CO 

H 

C 

H 

CO 

a 

<U 

oo 

M 

3 

<D 

00 

en 

•H 

CO 

fa 

CQ 

#> 

ON 

CO 

H 

■u 

CU 

oo 

29 


mmm 


S3 

w 

13 

h    o* 

w 


ecrccr 


c 

•H 

ca 

CD 

en 

CO 

PQ 

#. 

cn 

N 

=s= 

33 

S 

w 

.J 

CN 

P4 

• 

S 

OO 

< 

X! 

cr 

W 

cu 

u 

O 

Pn 

z 

H 

CD 

H 

en 

Pn 

cd 

H 

« 

PS 

W 

*« 

0 

•0- 

1 

ON 

CM 

CN 

H 

e 

ca 

CD 

0 

K 

CA 

3 

GO 

CD 

•H 

cn 

Pn 

m 

CQ 

«* 

vO 

CN 

rH 

4J 

CU 

W 

J 


30 


<r 


e.  <:  <  < 


« 


V 


v 


<^^W" 


> ,  ^ 


*T 


~7 


^O 


tt 


IT 


a 


w 
p 
cr 
w 
pi 
fa 


pa 


0) 

CO 

CO 

PQ 

<• 

ft 

^rfc 

N 

w 

w 

£ 

►J 

fa 

<r 

S 

rH 

<! 

X 

cr 

w 

cu 

rJ 

o 

fa 

ss 

H 

cu 

H 

w 

fa 

CO 

H 

pq 

HI 

w 

A 

O 

^O 

1 

00 

CN 

<n 

rH 

c 

CO 

O 

cj 

U 

en 

3 

60 

cu 

•H 

en 

fa 

CO 

pa 

A 

r- 

sO 

■u 

0) 

(/) 

31 


One  final  category  for  possible  signal  anomalies 
consists  of  five  different  phenomena  eacti  of  whica  occured, 
with  one  exception,  only  once  in  the  thirty-four  data  sets. 
Figure  14  is  an  example  wherein  three  different  "anomalies" 
were  present.  (It  should  be  noted  that  base  scan  = 
"various"  for  Figs.  14  and  15  means  that  each  of  the  three 
sections  of  the  waveforms  in  the  two  figures  hal  different 
base  scans  and  implies  that  this  in  no  way  detracts  from  the 
usefulness  of  the  figures  as  examples  of  the  affects  in 
question.)  Bins  1  througa  40  and  83  through  98  exhibit  a 
very  high  degree  of  signal  attenuation  not  seen  in  normal 
data  sets  (i.e.  data  sets  not  found  to  have  either  a 
receiver  malfunction  or  test  signal  present) .  Bins  41 
tnrough  87  exhibit  normal  signal  characteristics.  Bins  101 
through  200  evidence  a  combination  of  unusual  attenuation 
and  normal  signal  activity,  while  bins  231  througi  270  show 
either  a  general  distortion  of  the  signals  present  or  a 
fairly  high  amplitude  level  of  ambient  noise  and/or 
distorted  signals  plus  noise.  None  of  these  phenomena  could 
be  explained  by  receiver  aalf unction,  data  card  sequencing 
or  special  signal  characteristics  (i.e.  test  signal,  etc.). 
Figure  15  shows  a  data  set  exhibiting  no  continuity  of 
signal  whatsoever.  None  of  the  signal  "bumps"  can  be 
followed  in  a  uniform  manner  as  time  progresses.  The  theory 
that  a  serious  rsceiver  system  malfunction  was  tie  cause  of 
the  distortion  *as  borne  out  by  a  call  to  the  cognizant 
authority.  This  malfunction  was  also  verified  for  data  sets 
25  and  54.  The  sawtooth  affect  seen  at  the  aidlle,  right 
hand  side  of  Fig.  16  is  a  pattern  not  seen  elsewaere  either 
in  this  set  or  any  of  the  other  thirty-three  sets  analyzed. 
It  shows  wild  fluctuations  in  seven  adjacent  amplitude 
points  (2.8  kHz  apart)  ia  one  scan,  whereas  the  rsst  of  the 
scan  appears  to  be  normal.  No  explanation  for  this 
phenomenon  can  be  offerel.  Another  unexpected  anomaly  was 
the  two  flat-topped  signals,  centered  at  bins  124  and  128, 
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Figure  16  -   SAWTOOTH  PATTERN  EXAMPLE 
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found  in  Fig.  17.  It  appears  to  be  the  result  of  soma 
special  attenuation  applisi  to  those  frequency  bins  alone, 
because,  as  can  be  seen  by  the  rather  "large"  signal 
cantered  on  bin  178,  other  signals  appear  not  to  be 
affected.  The  additional  attenuation  applied  to  this  set 
during  reception  ,  10  d3,  is  lower  thin  levels  applied  to 
other  sets  which  were  not  similarly  affected.  The  best 
phenomenon  observed  in  this  final  category  is  f  ouad  in  Fig. 
18.  For  some  reason  (other  than  the  possibility  of  a 
different  form  of  receiver  malfunction  than  seen  before)  an 
otherwise  normal  data  set  lost  total  signal  continuity 
across  tne  entire  frequency  bandwidth  witi  subseqient  scans 
exhibiting  a  highly  distorted  nature.  While  tais  could  not 
have  been  caused  oy  a  card  sequencing  problem,  its  source 
cannot  be  identified. 


C.   ANOMALY  B  ANOMALIES 


The  Anomaly  B  progran  provided  a  different  perspective 
from  which  to  observe  tie  data  sets  since  it  presented  two 
dissimilar  plots  on  the  screen  at  the  same  time.  The  lower 
traces  consisted  of  five  Anomaly  k  style  three-dimensional 
plots,  this  time  only  eigaty  points  (or  bins)  wide,  while 
the  upper  five  traces  were  the  two-dimensional  time  vs. 
amplitude  plots  of  bins  seLected  from  the  traces  oelow.  It 
was  felt  that  this  different  perspective  would  enhance  the 
anomaly  identification  effort,  that  it  did  not  was 
disappointing. 

Of  the  eighteen  data  sets  analyzed  (3,  7,  16,  18,  20, 
21,  33,       36,    31,    38,  40,  41,  51,  53,  51,  80,  125,  and  139), 

Anomaly  B  proved  useful  in  only  two  instances:  (1)  to 
provide  a  graphic  presentation  of  the  three  signal  types 
used  in  the  Statistical  Aaalysis  Section  of  this  paper;  and 
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Figure  17  -   ATTENUATED  SIGNALS  EXAMPLE 
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(2)  to  illustrate  one  vary  interesting  phenomenon  which 
could  not  be  seen  by  the  Anomaly  A  presentation.  As  shown 
in  Fig.  19  this  phenomenon  consists  of  an  undulation  in  the 
amplitude  of  the  five  uppec  traces.  This  "lazy  sinusoidal" 
variation  appears  to  affect  some  signals  more  often  than 
others.  The  term  "lazy  sinusoidal"  should  be  usei  with  care 
since  it  is  not  meant  to  imply  tnat  there  is  any  strict 
periodicity  to  be  found  in  the  undulation  or  that  the  traces 
look  exactly  like  sinusoids.  However,  many  of  the  signals 
did  evidence  the  more  smooth  (and,  hence,  sinusoidal) 
variation  seen  in  trace  one  in  the  figure. 

Although  only  relatively  smooth  waveforms  ace  depicted 
in  the  figure,  this  effect  was  identified  over  a  large  range 
of  signal  types  including  those  whica  teyed  off  and  on; 
however,  it  was  much  easier  to  see  this  phenomenon  on  the 
smooth  signals.  This  variation  was  seen  in  every  data  set 
analyzed  although  it  occured  more  freguently  and  noticeably 
in  some  than  in  others. 

The  cause  of  tnis  undulation  could  be:  (1)  variations  in 
transmitter  output  power  levels;  (2)  receiver  system 
processing;  or  (3)  ionosphece-celat ed  phenomena  including 
Faraday  Rotation,  Sporadic-fi  and  fading. 
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Figure  19  -   AMPLITUDE  UNDULATION  EXAMPLE 
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IV.   STATISTICAL  ANALYSIS 


A.   DISCUSSION 


In  conjunction  with  the  graphical  analysis  approach  to 
identify  signal  anomalies  attrioutable  to  ionospheric 
properties,  an  attempt  was  made  to  glean  additional 
information  from  the  BSIGHAtf  data  by  means  of  statistical 
algorithms;  however,  it  should  be  noted  taat  no  attempt  was 
made  to  identify  anomalies  by  this  approach.  Bather,  tha 
programs  MIN/&AX  1  and  30/HAX  2  wars  developed  ia  an  effort 
to  determine  how  often  and  by  what  amount  (amplitude  in 
decibels)  both  signal  wavefronts  (scans)  and  individual 
signals  (bins  or  sets  d£  bins)  change.  It  was  hoped  that 
this  new  data  might  prove  useful  in  the  ievelopement  of  an 
algorithm  to  adjust  the  rnyciad  sensitivity  controls  found  on 
the  swept-tuned  receiving  equipment  used  to  acquire  and 
record  the  3RIGHAM  data.  Although  this  effort  aas  nothing 
directly  to  do  with  the  signal  analysis  effort  reported 
previously,  it  did  take  considerable  time  and  effort  to 
write  the  programs  and  create  the  desired  data  base.  In 
addition,  this  evolutioi  proved  extremely  helpful  in  tne 
overall  effort  of  understanding  signal  characteristics,  the 
primary  goal  of  this  stuly. 

For  the  purposes  of  tiis  inquiry,  three  different  types 
of  signals  (one  that  stayed  "on"  through  most  of  the  signal 
duration,  one  that  stayed  "off"  through  most  of  the  signal 
duration  and  oae  that  appeared  to  key  on  and  off  throughout 
most  of  the  signal  duratiDa)   were   chosen   from   data   sets 
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which  were  collected  it  various  tiies  of  the  day.  The 
quotation  aarks  are  used  (bare  and  throughout  tha  rest  of 
the  paper)  to  imply  a  ganeral  trend  notad  ovac  the  life  of 
the  signal  not  soma  unvarying  signal  characteristic.  Figure 
20  is  an  Anomaly  B  graph  plot  presenting  classic  sxamples  of 
the  three  signal  types.  Tha  "ON"  signal  (exhibiting  that 
"lazy  sinusoidal"  amplitude  variation  mentionai  in  tha 
graphical  analysis  sectioi)  is  traca  1,  the  "ON/DFF"  signal 
is  trace  2  and  the  "OFF'1  signal  is  trace  3.  Tracas  4  and  5 
are  merely  copies  of  traca  3  to  fill  out  the  plot.  A  table 
giving  all  partinent  data  on  tha  twenty-foir  signals 
analyzed  may  be  found  in  Fig.  21.  In  tha  figure,  sat 
frequency  refers  to  tha  base  frequency  of  tha  sat  when  the 
recording  was  initiatel  whila  101  maans  tha  time  of 
interception  for  the  data  sat.  kn  attempt  rfas  mada  to 
analyze  two  signals  froa  aach  of  four  diffarent  time 
categories  (interesting  from  the  ionosphere's  standpoint) 
, those  being  sunrise,  qodi,  sunset  and  aidnight.  With  one 
exception,  set  53  -  21D9LOCAL,  this  effort  was  realized. 


STATISTICAL  ANALYSIS  P&DGRAMS 


The  first  generation  algorithm  developed,  called  MIN/SAX 
1,  provided  for  a  scan-by-scan  (down  entire  frequency 
bandwidths  -  cross  sections  of  many  signals  -  at  fixed 
times)  analysis  of  a  variable  number  of  scans  witiin  a  givan 
set.  MIN/MAX  1  identifiad  eacti  of  the  peaks  and  valleys 
along  the  scans,  chose  the  highest  peaic  and  lowest  valley 
and  computed  the  average  peak  and  average  valley  heights. 
The  quantization  levels  given  oa  the  data  cards  ware  changed 
to  decibels  and  the  amount  of  attenuation  applied  at  the 
time  of  acquisition  was  added  back  in  to  gi/e  received 
signal  levels  accurately.  MIN/1AX  1  proved  to  b=  of  only 
small  value  since  signal  rfavefronts  vice  signals  themselvas 
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Figure  20  -   THREE  MIN/MAX  SIGNAL  TYPES 
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were  analyzed.   A  program  Listing  far  MIN/MAX  1  may  be  found 
at  the  end  of  this  paper. 

The  next  generation  effort,  tabbed  MIN/MAX  2,  focused  on 
specific  signals  (the  three  varieties  delineated  in  the 
discussion  section  above)  and  aided  a  number  of  more 
valuable  statistical  features.  Like  its  pc edecessoc , 
MIN/MAX  2  identified  each  of  the  peaks  and  valleys,  chose 
the  highest  peak  and  lowest  valley  and  computed  the  average 
peak  and  average  valley  heights. 

Its  next  feature  was  the  dateraiaation  of  twD  different 
types  of  signal  changes:  (1)  the  comparison  of  sample  point 
amplitudes  with  a  sliliag  refernce  point  to  determine 
long-term  changes  in  signal  levels;  and  (2)  tha  comparison 
of  adjacent  sample  point  amplitudes  to  determine  short-term 
changes.  The  former  was  initiated  with  the  bin  1  point  as 
the  reference  poiat  and  all  suoseguant  points  compared  with 
it  until  an  amplitude  was  found  which  differed  from  tne 
reference  by  a  previously  agreed  upon  number  Df  decibels 
(for  this  test,  4,  6,  12  and  18  dB  chaages  ware  used  in 
succession) .  That  point  then  became  the  reference  point  and 
the  process  continued  for  the  entire  2.4  minute  length  of 
the  signal.  The  number  of  sucn  long-term  changes  and  the 
elapsed  time  between  each  change  were  recorded.  Phis  data 
could  then  be  used  to  attain  the  goal  of  this  effort,  the 
gauging  of  how  often  and  oy  what  amount  the  various  signal 
types  changed. 

A  further  refinement  of  the  procedure  was  the  creation 
of  amplitude  histograms  of  the  elapsed  times  between  the 
changes  to  illustrate  the  effect  tnat  increasing  the  level 
of  decibel  fluctuation  which  constituted  a  change  had  on 
both  the  frequency  of  changes  and  the  distribution  of 
elapsed  times.  Since  the  receiver  sampling  rate  was  25  Hz, 
the  smallest  elapsed  time  possible  between  X  dB  fluctuations 
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was  0.04  seconds.  A  bcief  analysis  of  the  first  set  of 
elapsed  times  lad  to  tha  division  of  tha  time  line 
(abscissa)  into  twenty-one  time  bins  with  tha  inclusive 
boundaries  given  in  Fig.  22  and  labalad  N1  through  N21.  As 
a  further  refinement,  oins  N1,  N2  and  tha  first  segment  of 
N3  were  divided  into  individual  0.04  sacond  bias  labeled 
N1-1  through  N1-7,  N2-1  through  N2-3  and  N3-1  respectively. 
The  two  histograms  gave  aa  excellent  perspective  from  which 
to  watch  the  changes  in  signal  characteristic  as  the  amount 
of  decibel  fluctuation  tasted  was  stepped  from  '4  to  6  to  12 
and,  finally,  to  18  dB. 

The  final  output  of  MIS/MAX  2  was,  in  additisn  to  the 
peak  and  valley  statistics  mentioned  aoove,  a  listing  of  tne 
number  of  long-term  changes  that  took  place  in  =i:i  of  the 
thirty-seven  time  bins  of  interest,  the  percentage  of  the 
number  of  changes  under  each  time  bin  as  compared  with  the 
total  number  of  changes  (this  allowed  for  a  more  valid 
comparison  of  the  over-alL  change  picture  when  ai  increase 
in  decibel  deviation  ceguired  resulted  in  fe*er  changes 
tabulated  in  the  bins) ,  the  total  number  of  long-term 
changes  (called  NOM  in  tae  output)  and  the  total  number  of 
changes  in  adjacent  sample  points  (the  short-tera  changes 
mentioned  above).  In  addition,  the  mean  and  standard 
deviation  for  tae  elapsed  times  between  long-teca  changes 
were  computed.  A  listing  for  MIN/SAX  2  may  De  fDund  at  the 
end  of  this  paper. 


C.   OBSERVATIONS  AND  CONCLUSIONS 


There  were  a  number  of  surprises  which  resulted  from  the 
initial  computer  output  ccom  ttlN/MAX  2.  First,  the  high 
numoer  of  both  long-  and  short-term  changes  identified  was 
not  expected.   Out  of  the  possible  3597  changes  (there  ware 
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at  most  that  many  scans,  dc  sample  points,  available  for 
each  signal  of  interest)  the  number  of  changes  noted  in  some 
cases  reached  as  high  as  2000  for  the  lower  decibel 
fluctuation  levels.  Second,  the  elapsed  time  histograms 
tabulated  showed  a  marked  tendency  to  be  heavily  damped 
exponentials  in  nature  with  the  0.0'4  second  bin  (N1-1) 
having  the  highest  number  of  occurences  throughout.  And 
third,  there  was  a  high  degree  of  overlap  between  the 
statistical  characteristics  generated  foe  the  three  signal 
types  chosen  indicating  that  there  was  no  absoLute  way  to 
categorize  a  signal  as  "OS",  "ON/OFF"  or  "OFF"  merely  by 
looking  at  its  statistics.  What  was  expected  (and  lid 
occur)  was  that,  as  tha  number  of  decibels  used  for 
comparison  increased,  the  exponential  distributions  "flowed" 
to  the  right  indicating  gceater  elapsed  times  between  larger 
decibel  fluctuations. 

Figures  23  -  29  represent  the  bulk  of  the  meaningful 
data  garnered  from  the  MIN/SAX  2  application  to  the 
twenty-four  signals  of  interest.  In  speaking  of  the 
characteristics  which  differentiate  the  three  arbitrarily 
chosen  and  defined  signal  types,  Figs.  23  -  28  provide  those 
statistics  which  can  be  of  value  in  coming  to  a  logical 
conclusion  concerning  tneir  differences.  As  would  be 
expected  from  "ON/OFF"  signals,  they  showed  a  greater 
difference  between  both  average  peak  and  valley  values  and 
between  "maxmax"  and  "minnin"  values  than  did  the  other  two, 
implying  more  signal  level  fluctuation.  The  "OFF"  signals 
consistently  snowed  ths  smallest  fluctuations  as  would  be 
expected  from  a  series  of  either  highly  decayed  signals  or 
high  but  fluctuating  ambient  noise.'  Based  on  these 
statistics  alone  it  woald  be  nearly  impossible  to 
distinguish  with  a  high  degree  of  certainty  between  an  "ON" 
and  an  "ON/OFF"  signal  given  only  the  parameters  in  those 
figures.  Figure  29  provides  both  a  summary  of  the  data 
concerning  tne  number  of  ohanges  tabulated  for  each  decibel 
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KEY:     PEAK  =  AVERAGE  OF  ALL  MAX  PTS  IN  EACH  BIN  (SIGNAL) 

VALLEY  =  AVERAGE  OF  ALL  MIN  PTS  IN  EACH  EIN  (SIGNAL) 

NAXMAX  =  LARGEST  PEAK  VALUE 

MINMN  =  SMALLEST  VALLEY  VALUE 

PVR  =  RANGE  OF  CIFFERENCES  BETWEEN  PEAK  £  VALLEY 

PVA  =  AVERAGE  CIFFERENCE  BETWEEN  PEAK  £  VALLEY 

K4R  =  RANGE  OF  CIFFERENCES  BETWEEN  MAXMA*  6  MINMIN 

MMA  -  AVERAGE  DIFFERENCE  BETWEEN  MAXMAX  £  MINMIN 

PEAK  AVG  =  AVERAGE  CF  ALL  PEAKS  FCR  "CN/OFF"  SIGNALS 

VALLEY  AVG  =  AVERAGE  CF  ALL  VALLEYS  FCR  "ON/CFF"  SIGNALS 


Figure  24  -   "ON/OFF"  SIGNAL  PEAK-VALLEY  STATISTICS 
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KEY:     PEAK  =  AVERAGE  CF  ALL  MAX  PTS  IN  EACH  BIN  (SIGNAL) 
VALLEY  =  AVERAGE  CF  ALL  MIN  PTS  IN  EACH  BIN  (SIGNAL) 
NAXMAX  =  LARGEST  PEAK  VALUE 
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PVR  =  RANGE  OF  CIFFERENCES  BETWEEN  PEAK  6  VALLEY 
PVA  =  AVERAGE  DIFFERENCE  EETWEEN  PEAK  £  VALLEY 
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l*MA  =  AVERAGE  C  IFFERENCE  .BETWEEN  MAXMAX  £  MINMIN 
PEAK  AVG  =  AVERAGE  OF  ALL  PEAKS  FCR  "OFF"  SIGNALS 
VALLEY  AVG  =  AVERAGE  CF  ALL  VALLEYS  FCR  "OFF"  SIGNALS 


Figure  25  -   "OFF"  SIGNAL  PEAK-VALLEY  STATISTICS 
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KEY*:   D6   =  D8  CHANGE  TESTED 

NUM  =  NUMBER  OF  LONG-TERM  CHANGES 

ADJ  =  NUMBER  OF  SHORT-TERM  (AOJACENT  POINT)  CHANGES 

N_   =  PERCENTAGE  OF  LONG-TERM  CHANGES  IS  TIME  BIN  N. 


Figure  26  -   "ON"  SIGNAL  ELAPSED  TIME  STATISTICS 
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KEY:   CB  =  CB  CHANGE  TESTEC 

NLM  =  NUMBER  CF  LONG-TERM  CHANGES 

ACJ  =  NUMBER  OF  SHORT-TERM  (ACJACENT  POINT)  CHANGES 

N_  =  PERCENTAGE  CF  LCNG-TERN  CHANGES  IN  TINE  BIN  N 


Figure  27  -   "ON/OFF"  SIGNAL  ELAPSED  TIME  STATISTICS 
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KEY:   CE   =  CB  CHANGE  TESTED 

NUM  =  NUMBER  OF  LONG-TERM  CHANGES 

ACJ  =  NUMBER  OF  SHORT-TERM  (ADJACENT  POINT)  CHANGES 
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Figure  28  -   "OFF"  SIGNAL  ELAPSED  TIME  STATISTICS 
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Figure  29  -   OVER-ALL  dB  VS.  CHANGE  DATA 
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comparison  level  and  a  general  means  of  distinguishing 
between  the  three  types  of  signals.  Althougn  there  are 
definite  differences  between  the  average  figures  for  the 
various  dB  levels  and  some  differences  in  the  ranges  of 
values,  the  overlap  encountered  in  the  ranges  mates  signal 
differentiation  by  this  means  risky.  In  snort,  what 
distinguished  the  "ON"  fcoii  the  "ON/OFF"  and  "3FF"  signals 
was  thoroughly  subjecti/e  and  only  loosely  corroborated  by 
the  statistics  generated. 

Additional  information  can  be  retreived  from  Figs.  26  - 
28  which  provide  the  number  of  changes  (both  Long-  and 
short-term)  found  when  the  iecibel  differential  wis  "walked" 
from  4  to  18  dB  and  the  percentage  of  long-teem  changes 
which  fell  into  the  most  meaningful  of  the  tairty-seven 
histogram  bins  (for  elapsed  times  between  long-term 
changes)  .  As  is  clearly  shown,  the  number  of  caanges,  NtJM 
and  ADJ,  decreased  as  the  dB  level  increased  while  the 
percentage  figures  in  columns  N1  and  N1-1  evidenced  a 
movement  away  from  the  first  elapsed  time  bin.  Both  of 
these  characteristics  imply  that  the  average  time  between 
signal  level  changes  at  the  various  do  settings  increases  as 
the  number  of  dB  increase.  To  be  specific,  by  taking  the 
data  from  Fig.  30  it  could  be  calculated  that  the  average 
time  elapsed  between  4  and  6  dB  f luctuat uions  foe  the  three 
signal  types  were: 
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Figure  30  -   CURVE-FITTING  DATA 
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This  data  confirms  tiat  the  "ON/OFF"  signals  fluctuated 
more  rapidly  than  the  other  two  categories  while  the  "ON" 
signals  proved  to  be  more  stable  than  the  "OFF". 

These  results  imply  that,  if  the  agency  that  operates 
the  swept-tuned  receiving  system  used  to  acquire  and  record 
the  BRIGHAM  data  desires  to  acquire  FSK-ty?e  signals 
(roughly,  the  "ON"  variety),  it  could  get  away  with  changing 
the  sensitivity  adjustments  much  Less  frequently  than  if 
Morse-type  signals  (rougaly,  the  "ON/OFF"  variety)  were 
desired.  This,  plus  the  rough  data  provided  in  Fig.  23, 
could  provide  a  good  data  base  from  which  a  dial  controlling 
algorithm  could  be  adjusted  to  meet  the  needs  of  the  signal 
types  of  interest  at  any  tine. 

As  was  mentioned  above,  MIN/MAX  2  provided  for  the 
calculation  of  the  mean  and  standard  deviation  of  the 
elapsed  times  tabulated  foe  eaca  of  forty-eight  signals  (the 
twenty-four  signals  of  interest  at  both  '4  and  6  dB 
comparison  levels).  Although  the  aean  values  were  accurate 
the  standard  deviation  values  were  not  since  tae  formula 
used  in  their  calculation  is  valid  for  only  normal  (or 
Saussian)  distrib Jtions.  Since  the  elapsed  time  histograms 
were  apparently  exponential  in  nature,  another  form  of 
variance  determination  was  required.  This  effort  called  for 
a  curve-fitting  technique  to  identify  the  histograms  as 
being  of  a  particular  distribution  type,  use  of  the 
chi-squared  test  for  gooliess-of -f it  and  manipulation  of  the 
chosen  distribution's  raeaa  and  variance  formulas  to  derive 
the  information  desirei.  Attempts  were  made  to  fit  three 
distribution  functions  to  the  histograms  (3amma,  Poisson  and 
Oeometric)  with  the  3anma  function  proving  to  be  the  most 
successful.  The  Gamma  distribution  function,  with 
parameters  C<,/^,  0<(X  <<=**,  0</3   <°°  is  defined  as  follows: 
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f(x)  = 


0, 


x<0 


Where,  Mean  =  <*/3  ,  Variance  =c*/3  . 

An  algorithm  was  written  which  read  in  tha  observed 
histogram  frequencies,  calculated  the  expected  Gamma 
frequencies  (using  an  International  Mathematical  Subroutine 
Library  subroutine  called  MDGM1)  and  calculated  the 
chi-sguared  value  (essentially,  the  sum  of  the  least  squared 
differences  betwsen  the  observed  and  expected  fceguencies) 
for  use  in  the  goodness-of-f it  test.  The  chi-sguired  value 
was  then  compared  wita  the  appropriate  value  found  in  the 
chi-sguared  distribution  table  (appropriate  as  to  degrees  of 
freedom  and  significance  level);  if  it  fell  below  that  table 
entry,  the  fit  was  considered  good,  if  not  the  fit  was  bad. 
As  discussed  in  Refs.  7  and  8,  a  good  rale  of  thamb  is  that 
the  smaller  the  chi-sguared  value  for  a  givea  number  of 
degrees  of  freedom,  the  better  the  fit.  ks  a  word  of 
caution,  even  an  excelleat  fit  (extremely  small  rhi-sguared 
value)  does  not  guarantee  that  the  distribution  is  really 
the  one  oeing  tested,  it  narely  states  that  the  data  do  not 
present  sufficient  evidence  to  contradict  the  hypothesis 
that  the  histograms  possess  Gamma  distributions. 

For  the  curve- fitting  at  hand,  knowledge  of  tie  mean  and 
an  iterative  variation  in  c<  peraitted  an  iterative 
calculation  of  /5  thus  providing  all  of  the  necessary 
information  for  tie  algocitnm  discussed  above.  Tie  process 
was  refined  until  the  smallest  chi-sguared  valae  could  be 
identified.  It  saould  be  noted  that  only  the  N1  through  N15 
time   bins   were  used  since  very  few  elapsed  tines  ever  fall 
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beyond  five  seconds  and  that,  as  tha  individual  histogram 
composition  dictated,  either  seven  or  three  bins  were 
actually  used  in  the  algorithm  (following  a  suggestion  that 
bins  containing  five  or  fewer  counts  be  combined  to  enhance 
the  accuracy  of  the  test) .  The  bin  limits  impliei  that  the 
degrees  of  freedom  used  were  five  oc  ona  respectively. 

Of  the  forty-eight  histograms  to  be  fitted,  six  had 
counts  in  too  fen  bins  to  be  useful  (would  have  inplied  zero 
degrees  of  freedom  -  not  found  in  the  taole) ,  eigateen  were 
bad  fits  and  tha  remaining  twenty-four  provided  good  fits  to 
the  Gamma  distribution.  Figure  30  gives  the  rasults  of  the 
curve-f itting  effort.  Knowledge  of  <2<  and  tha  lam  implias 
knowledge  of  {3  which,  ii  turn,  implies  (by  definition) 
knowledge  of  the  vaciaice.  K  combination  of  tie  mean  and 
variance  values  should  be  useful  in  the  determination  of  how 
often  the  sensitivity  dials  on  tha  receiving  system  should 
be  adjusted  for  the  dasirad  signals  to  be  acquired. 

Figure  31  provides  a  breakdown  of  the  c<  values  vs. 
signal  type  and  amplitade  change  variable  for  tha 
twenty-four  histograms  which  "fit"  the  Samma  distribution. 
It  was  felt  that  tais  would  permit  the  grouping  of 
histograms  with  raspect  to  some  combination  of  the  three 
categories  above.  The  six  ranges  shown  in  the  figure  imply 
that  any  grouping  by  q<.  alone  would  be  tenuous  indeed  as 
they  overlap  quite  heavil/.  The  only  discernaQle  trend  is 
the  apparent  lowering  of  the  average  cK.  figare  as  the 
amplitude  change  variable  increases  from  4  to  5  dB.  This 
would  seem  to  contradiot  the  defining  equations  for  the 
3amma  distribution  whioa  inply  that  a  lower  c<  implies  a 
more  "damped-exponential''  form  for  the  distribution  (when, 
in  fact,  the  more  damped  case  is  found  at  the  4  d3  change 
variaole) .  However,  tha  shape  of  the  distributioa  is  also  a 
function  of  /3  which,  as  defined  in  Sean  =  CK/5  ,  tends  to 
affect  the  distribution  saape  in  a  aore  radical  fashion 
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4  DB 


«CN"  "ON/OFF"  "CFF" 

SE7/eiN    ALPHA      SET/BIN    ALPHA      SET/BIN    ALPHA 

8C/82     0.18        21/273    0.85         3/169  0.87 

126/57     C.72  7/91  0.34 

20/259  0.94 

53/5C  1.05 

80/73  1.05 

RANCE:   0.18-0.72    RANGE:   NONE         RANGE:   C. 34-1. 05 
MEAN:   C.45  MEAN:   0.85  MEAN:   C.85 


6  DB 


«CN"  "ON/OFF"  "OFF" 

SEIZE  IN    ALPHA      SET/BIN    ALPHA      SET/BIN    ALPHA 

7/8      0.18         3/154    0.93         3/169  0.34 

2C/265    C.16        21/273    0.49         7/91  0.18 

21/15C    0.045       53/8      0.065       20/259  C.50 

126/5 1     0.62        8C/60     0.18        21/187  C.44 

126/74     0.28        53/5C  C.28 

81/11  C.76 

126/8  C.51 

RANGE:   0.045-0.62   RANGE:   0.065-0.93   RANGE:   C. 18-0. 76 
MEAN:   C.25  MEAN:   0.39  MEAN:   C.43 


Figure  31  -   SIGNAL  TYPE  VS.  ALPHA  TABLE 
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than  does  (X  .  in  fact,  as  tae  change  variable  increases, 
the  mean  increases  (see  Fig.  30)  ana  <X  seems  to  decrease 
(see  Fig.  31)  implying  a  very  pronounced  increase  in  ^ 
which,  in  the  iefining  sguations,  is  found  to  be  most 
prominent  in  the  exponential  teem;  but,  this  implies  a 
distribution  shape  less  heavily  dampled  in  the  lorfsr  (4  dB) 
change  variable  level.  Tiis  trend  of  decreasing  values, 
then,  was  to  be  expects!  and  appears  to  be  the  only 
meaningful  product  of  the  analysis  of  the  Fig.  31  data. 


As  an  aside  on  the  snoject  of  signal  continul 
the  brief  recoried  time  span  of  the  waveforms  ( 
minutes) ,  it  should  be  noted  that  data  set  number 
("ON/OFF")  ani  9t  ("OFF")  exhibited  vary 
characteristics  when  compared  with  the  other  s 
their  classes.  As  sno*n  on  Fig.  27,  bin  81  stil 
1000  changes  at  13  IB  iiffecence  while  the  other 
signals  had  fallen  off  to  much  smaller  numbers, 
that  bin  had  two  changes  at  54  dB  when  ever- 
numbers  were  applied.  rhs  "OFF"  signal,  bin  91,  3 
similar  degree  of  consistency  as  it  had  many  fa 
changes  than  the  other  "OFF"  signals  wnile  fallii 
more  quickly  as  the  number  of  decibels  increased 
28) .  In  fact,  this  signal  had  only  a  few  more  t 
dB  changes  implying  exceptional  amplitude  contin 
compared  with  any  of  the  twenty-three  oths 
analyzed.  To  give  a  mora  grapnio  representation 
phenomenon,  Figs.  32  and  33  compare  the  long-t 
changes  found  in  data  sets  3  and  7  with  the  latts 
to  be  the  more  consistant  of  the  two,  although 
nas  a  good  measure  of  uniformity  throughout  its  36 
The  base  scans  for  each  of  the  five  soan  sets  ( 
the  beginning,  niidle  and  end  of  each  data  set  - 
are  provided  at  the  rignt  of  the  pictaces. 
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Figure  32  -   LONG-TERM  SIGNAL  CONTINUITY  -  SET  3 
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Figure  33  -   LONG-TERM  SIGNAL  CONTINUITY  -  SET  7 
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V.   COMMENTS  AND  RECQNMEN DATI3NS 


Of  the  two  signal  anomaly  ideatif ication  techniques, 
Anomaly  A  proved  to  be  the  most  fruitful,  NarrDw  and  wide 
wavefront  discontinuities,  fading,  frequency  saifting  of 
entire  wavefronts,  frequency  shifting  of  individual  signals 
and  several  miscellaneous  varieties  were  seen.  rwo  of  the 
above  (discontinuities  aid  shifting)  i/ere  reported  on  in 
Ref.  1  with  the  rest  being  unique  to  this  study.  Several 
examples  of  pseudo-anomalies  were  included  to  perform  the 
dual  function  ot  warning  against  the  false  assimption  of 
anomalous  performance  siaply  because  the  plot  seemed  to 
imply  it  and  to  provide  pictures  with  the  narrati/e  in  such 
cases.  Analysis  with  Anomaly  B  proved  to  be  relatively 
unproductive  as  only  the  interesting  and  frequently  seen 
phenomenon  of  signal  amplitude  undulation  (reported  in  Ref. 
1  and  expanded  upon  here) =  oeing  noted.  The  diverse  and 
amazingly  complex  natnre  of  the  ionosphere  and  its 
interaction  with  its  environment  precludes  tie  precise 
assignment  of  a  cause  to  every  effect  (i.e.,  anomaly). 
Suffice  it  to  say  that,  except  for  those  anomalies  which  can 
be  explained  away  by  equipment  malfunction,  nothing  has  been 
found  in  this  effort  whicn  will  simplify  the  model  of  the 
ionosphere  as  an  ever- changing  Helium  in  waici  literally 
anything  can  happen  to  a  radio  signal  which  travecses  it. 

The  statistical  analysis  work;  performed,  from  creation 
of  the  MIN/MAX  algorithms  to  the  tedious  curve-fitting 
effort,  resulted  in  a  wealth  of  predominantly  raw  data 
wnich,  with  the  mean  and  variance  values,  should  provide  a 
foundation  upon  whici  a  receiver  systea  sensitivity 
dial-controlling  algoritna  could  be  refined.   It   was   noted 
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in  the  curve-fitting  section  that  it  is  perhaps  the  weakest 
link  in  the  statistical  study  since  the  chi-33 jared  test 
guarantees  neither  that  the  elapsed  time  aistograms  are 
3amma  distributions  nor  that  thera  night  exist  other 
distributions  which  woald  "fit"  the  histogram  data  as  well 
as  the  Gamma  function  did.  However,  a  statistically 
reasonable  argament  can  be  lade  that  the  histograms  (at 
least  half  of  them)  are  pcobaDly  Saitna  in  nature  and,  hence, 
the   variance  values  calcalated  can  ba  accepted  as  accurate. 

Should  additional  work  be  done  with  the  BRI^aAM  data, 
the  following  ramarks/r econ  nendations  could  prove  useful: 

(1)  all  of  the  complata  data  sats  in  hand  nave  been 
thoroughly  anaLyzed  witn  program  Anomaly  A  -  additional 
effort  along  this  line  mijat  be  a  waste  of  tima  and  energy; 

(2)  althouga  nearly  twenty  of  taa  data  sets  in  hand  have 
not  been  analyzed  with  program  Anomaly  3,  the  disappointing 
results  thusfar  might  act  as  a  warning  as  to  expected  future 
rewards  from  additional  affort; 

(3)  if  Anomaly  A  and/or  Anomaly  B  are  to  oe  ased  again, 
the  6AKSCN  featura  found  ia  tha  NErA5Y!lBDL  subroutine  (it 
provides  for  taa  backing  ovar  of  a  specified  number  of 
files,  or  data  satsj  shoaLd  be  aada  to  work  as  a  time/labor 
saving  device  -  the  problem  is  in  the  magnetic  tape  unit, 
not  the  program; 

(4)  the  recommendation  (aade  initially  in  Raf.  1)  for 
the  creation  of  a  compatsr  algorithm  to  searca  out  the 
anomalies  in  an  automatic  and,  henca,  more  rapid  nanner  than 
the  visual  approach  is  guardedly  endorsed  herein  -  although 
the  speed  at  waich  the  data  sets  could  be  analyzsd  would  ba 
improved  oy  at  least  an  order  of  magnitude,  tae  extreme 
complexity   of   trying  to  tell  the  computer  what  to  look  for 
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as  well  as  the  loss  of  that  very  intangible  but  valuable 
human  discretions ry  power  could  conbine  to  male  this  effort 
a  quagmire; 

(5)  should  additional  work  be  done  utilizing  the  MIN/MAX 
routines,  Anomaly  B  shouli  De  used  first  to  eihance  the 
accuracy  in  selecting  tae  different  signal  types  to  De 
analyzed  -  that  this  was  lot  done  in  the  present  study  may 
have  resulted  in  the  rather  large  ranges  found  in  the  data; 

(6)  the  accuracy  of  aiy  future  car ve-f itting  rfork.  couli 
be  enhanced  by  the  incLasion  of  the  tfeibull  and  Lognormal 
distributions  (or  any  of  the  level- crossing  I istr ibutions) 
to  the  list  of  possibLe  curves  to  which  the  eLapsed  time 
histograms  are  Batched  because  sets  of  such  lata  can  be 
fitted  to  more  than  one  lis tribation  with  equal  accuracy  and 
validity  -  this  was  not  lone  here  oecaase  the  suggested 
distrioutions  were  complex  eaoaga  and  the  tins  remaining 
short  enough  that  a  gooi  effort  couli  not  have  been  put 
forth. 


67 


GLOSSARY 


1 •  Sgoradic-E:  thin  horizontal  layer  (oc  patent  of  high 
electron  density  enbedded  in  the  regular  E  layer 
resulting  in  signal  deflection  and/or  absorption. 

2-  Nig^ht^Dax  Asy_iii[netry_:  Lass  attenuation  at  nigat  than  lay 
due  to  difference  in  the  ionizing  affect  o£  the  sun 
resulting  in  discontinuities  airing  the  sunrise/sunset 
periods . 

3.  Geomacrnetic  N25.zl§.2t2£22iSY:  different  attenuation 
effects  on  signals  travelling  wast  to  east  rather  than 
east  to  west  due  to  the  earth's  magnetic  fieLd  lines. 

**•  E.2liJL  Qa£  Absorption  X£3Aj_:  absorption  of  radio  signals 
passing  the  polar  regions  -  due  to  a  combination  of 
solar  flare  activity  and  the  shape  of  the  earth's 
magnetic  field  lines. 

5-  Au£2£§.l  DisDla^:  brigatly  colored  northern  latitude 
phenomena  disruptive  to  communications  -  also  as  the 
result  of  a  combination  of  solar  flare  activity  and  the 
earth's  magnatic  field  lines. 

6*  £araday_  Rotation:  any  Linearly  poLarized  wave  travelling 
in  the  direction  of  a  magnetic  filed  results  in  its  two 
circularly  polarizei  components  travelling  at  different 
velocities  and  thus  the  plane  of  polarization  will 
rotate  with  distanca. 
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7.  Ducting:  if  a  radio  wa/a  comas  into  contact  with  a  region 
of  inhomoganaoas  refractive  indices  the  wa/e  could  be 
trapped  between  t  *o  layers  and  guided,  as  in  a  leaky 
wave  guide,  away  fron  its  intended  iestinatiDu. 

8-  £§.iliD:3.  !D3lkiSli£  ^ffa_2ti_:  sadden  ionospheric 
disturbance  (S.  I.  3.)  which  producas  a  coapLate  radio 
"fade  out"  lasting  from  a  faw  minutes  to  an  hour  or 
more   -  caased  by  soLar  flare  activity. 
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aa  .  c 


flde  f i  ne 

XINT 

380 

<* de  f  i  ne 

YINT 

8000 

*de  f i  ne 

DELY 

3200 

Ode  f  i  ne 

DELX 

16 

tfde  f  i  ne 

AMPY 

8000 

tfde  f  i  ne 

NLINE 

1 

Kde'ine 

PI 

3.1«1596SaS 

#de  f  i  ne 

NPT 

150 

i  nt  zero 

o; 

i  n  t  one 


l; 


st  rue t  { 

i  nt  not • 

i  n t  line» 

i  nt  1  dx ; 

int  ldy,' 

i  nt  cot r  ; 

int  nsot; 

\ nt  soot  1  » 

\  nt  soot  2'f 

i  nt  spot  3; 

int  h  r  »  m  i  n  »  sec  » 

i  nt  n  f  t  ; 

1  nt  1 aq» 

i  nt  sr! 

int  Idl3)  ; 

i  nt  lo» 

int  set; 

int  mo  < day  t  y  r ; 
>   head»   *h ; 

int  hbuf  (3001  ; 

st  rue  t   data     < 

int  x ; 

int  yd; 

)  ; 

st  rue  t  dat  a  *d> 
int  dbuf  [40001  ', 
int  np; 
int  nip  8  ; 


int  idev#pdev,sclx»scly#hiasx,tdev; 

mai  n  ( arac  > argv ) 
int  *  *arqv ; 
1 

int  i  >  )  t n; 

char  *cs; 

if((tdev  =  open ( "/dev/spp" , 1) )  <  0)< 

print f ("cannot  ooen  spp  »"); 

exi  t  (  )  ; 
) 
if((pdev  =  open("/rtev/rvo"»l))<0){ 

pri  nt  f  (  "  cannot  open  rvn  tt"  )  ', 
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e  x  i  t  (  )  ; 

> 

if((idev  =  open( "/dev/rmt 6" # 0 ) ) <0 )  ( 

prfntf  ("cannot  open  rmt6  *>")', 

e«i t  (  )  ; 
) 

i  f (arqc  >  1  )  ( 

cs  =  arqv  f 1 )  ? 

n  =  o; 

while  (  *cs  >=  '0'  *&  *c<;  <=  '9') 

n  :  n  *  10  ♦  *es  t  +  -  '0'; 
n  =*  1  l ; 

for  (  i  =  0  ;  i<n;  i  ♦  +  ) 

inpUdev,  dbuf#f00); 

printf  ("number  of  records  skipped,  XdU",n)i 
) 
while((inp(idev,hbuf,B00))>0M 

h  =  hbuf; 

j  =  o; 

f  o  r  (  i  =  0 ;  i  <  h  ->  line;  i  +  +  )  ( 

n  =  inp( 1dev,fcdbuf  ( j 1  ,  BOO); 

j  =  ♦  2* (h->npt ) ; 

1 

sclx  =  030;  scly  =  03; 
bi  asx=50; 

np  :  h->npt  *  h->liner 
sea  1 e (  ) ; 

plot  (); 

c vers (pdev , 020  )  ; 
st  t  v (pdev ,  Jone ) ; 
1 


#de  f  i  ne 

NBYT   2hH 

#de  f i  ne 

NBLK   0 

Ode  H  ne 

NSL   1250 

tfde  f  i  ne 

DRAW       1 

st  rue  t  i  pt  < 

i  nt 

y ; 

char 

*xpb; 

char 

y  i  nc ; 

char 

xdi  r ; 

fnt 

yf ; 

i  nt 

cxp; 

Int 

Mllnk! 

i  nt 

*blini<; 

1     i  tab  [7*NBYTJ  ,  Mo,  ♦  i  a  J 


char    ob  (NBYT1  ; 
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int  *dp  12000)  ; 


plot  () 


//clear  plot  buffer 


//move  to  top  of  plot  area 
//plot  scan  line 


st  rue  t  i  pt  *s  ; 
int  i  t  j  t  i  e» 

sort  ().* 

for(i  =  0;  i  <  MBYT;  i+t) 

pom  =  o; 

for(i  =  0;  i  <  NBLK;  i+O 

wri  te(pdev,ob,MBYT ) ; 
ip  =  0;  ia  =  itab; 

j  =  o; 

fori  i  =  nsl;  i  >  -l :  «  — )  < 

while(j  <  nv    88  *dp(j]  ==  i)( 

ie  =  sip(dptj))'  //set  up  plot  point  for  interpolation 

i  f (ie  ==  1 ) 

return; 

nib(i)»  //set  up  plot  buffer 

wr i t e (nde v » ob t NR Y T ) ; //p 1 ot  line 

> 


int  i  sor t 


sort ( ) 
< 


regi  st  er  i  ,  k  ,  t  ', 


i=0i 


for(d  =  dbuf;  d  <  8dbuMnp«2J;  d*  +  ) 
dp ti ♦♦]   =  R(d  ->  yd); 

k=np  ; 

whi le     (     k    =>>     1     ) < 
J  sort**; 

while     (     j  sor t     ) ( 
jsort     =    0; 
for    (     i=0;     i<    Cnp-k) J     i  +  O 

i  f     (     *dD  I  i 1     <    *dp  (i  +k )     )  ( 
t=dofi]  ; 
dp  (i  1  =do  H+kl  ; 
dp  t  i  +  it  1  =  t  ; 
j  sort  ♦  +  ; 
> 


I 


1 
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s  ip(dp  i ) 

i  nt     *dpi  » 
< 

int     xi  ,  y  i  ,xl ,  y  I  , xr , yr ,  incyr fyrdi r»,cx(»s, lex,  inc! 

i  nt     \  t  i  e  J 

s  =  dpi  ; 

lex    =    1  ; 

yi     =     *dpi--?     xi     =     *dpi; 


1 f (  ♦  ♦  s    <    &(dbuf  (2»np)  )  )  < 
x  r    =    *  s  t  +  ; 
y r    = '  *s  ; 
i  f  (yr    8.    DRAW)  { 

if((inc     =    y  i     -    yr  )     >=    0)  < 
i  f  (  ( incy    =     inc )     ==    0) 

dl r«     =     nip; 
e  1  se  < 

f  e>  r  (  i     =     0;     ((incy     =     (  i  n  c  /  (  n  i  p    >>     i))) 
d  i  r  x     =     (I     <  <     i  )  ; 
) 
f  y    =    y  p  » 

ex   =   o?oo; 
lex   =   o; 

ie    =     stacMvi-incy,incy»fy,xi,dirx,cx); 
i  f (i     ==     1  ) 

re  f  urn ( 1  ) ; 


=  =    0)  ;     i ♦  +  )  ; 


} 


> 


> 

i  f ( — dpi     >=    dbuf ) < 
i  f (yi     &    DRAW)  < 

y I     =     *dp  i -- ; 
x 1     =    *dp  i ; 
i  f  (  ( i  nc     5    y  i 
i  f ((incy 


yl  )  >=  0)  < 
inc)  ==  0) 
-nip» 


d  i  r  x 
e  I  se  { 

for(i     =0;     ((incy    =     (inc/(nip    >>     i)H     = 

dirx     =    -(1     <<     i); 
) 

f  v    =    y  I  ; 
ex    =    I c  x ; 

ie    =    s t ac k ( y i - i ncy ,  i ncy t  f y ,  x i -  1 , di r x , c x ) ; 
i  f  (  i  e    =  =     1  ) 

return! 1 ); 


=    0  ) ;    i  +  +  )  ; 


> 

return(O)  ; 


stack(arbfC>ddfef  f) 

int     a  ,\n  r  c  ,  <lci,p  ,  i  ', 
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{ 


*s; 


1nt 

int  i  ; 

st  rue t  i  ot 


2  i 


i  a  -  >  y  =  a ; 

i  a  ->  vine  =  b ; 

i  a  ->  y f  =  c; 

i  f (dd  >=  NBYT  li    dd  < 

pr i  nt  f ( "bad  scale 

exi  t  ( ) ; 
} 

i  a  ->  xpb  -    dd  ♦  obi 
ia  ->  xd  i  r  =  e ; 
ia  ->  cxp  =  f; 


0)( 

x=%dS",dd) ; 


i  f (ip    ==    0) { 
ip=  i  tab; 

ip->fHnk  =  ip->bHnk  =  Oi 
i  a  ♦  ♦  ; 

ia->bl i  nk  =  i  p; 
i  a->f  Hnk  =  0; 
> 
e  I  se  < 

s=1a->bl ink; 
s  -  >  f  1 i  n  k  =  i  a  ; 
i  f  (  i  a->f 1  ink  ~-     0)  < 
s  =  iai 
j  f ( ♦+  i  a  >r  Ri  tab(7*NBYT]  ) < 

print  f ("i  tab  overflow  »") J 
ret urn ( 1 ) J 
> 

ia->b1 i  nk  =  s ; 
i  a->f 1  ink=0; 
> 
e  i  se  { 

s  =  ia  ->  flink? 
s  ->  blink  =  ia? 
ia  ->  flink  =  0; 
i  a  =  s  » 


> 


> 


nib(sl ) 
int 


si  ; 


int  i  f  j i n ; 

int  *S ;  S  =  ioi 

w  h  i  1  e  (  s  )  ( 

i  f ( s  ->  y  f  <0  ) 

i  f ( (s=f ree(s)  )  ==  0) 
return; 
*s  ->  xpb  =s  ->  cxp; 
i  f ( s 1  ==  s  ->  v  )  < 

\      :  ((n  :  s  ->  «rtir)  >0 


//setuo  plottinq  buffer 
//delete  point 


// 


i  ne  break 
-n)  ; 
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i  f  (  n  <  0  )  //left 

f  o  r  ( i  =  0 ;  j  <  i ;  J  t  + )  < 
i  f  (s  ->  c*r>  ==  0)< 
*s  ->  xnb  =♦  I  »* 
s  ->  c  xp  =  1  ; 
> 

»s  ->  xnb  =!  (s  ->  cxp  =<<  1); 
> 
el  se 

f  o  r  (  j  =  0  ;  j  <  i  ;  j  *  f  )  { 

*S->XOb  =!   (s->Cxp  =>>   l)r 

s->y  =-  s->yinC'* 


> 

i  f (si  <=  s->vf) 

s->y f  =  -  1  » 
s  =  3*>n  ink! 


/ /end  o  f  po i  n  t 


f  ree ( s ) 

int  *s  J 
< 

int  *t ; 
int  1  ; 
st  rue  t  i  pt  *  z  ; 

*s  ->  xpb  =  o; 

i  f (s->b!  ink  ==  0)  < 
ip  =  s->  f link i 
ip->bl ink  =  0; 

t  =  i  p ; 

> 

e  1  se  { 

t  =  s->Hl ink; 

t ->f 1  i  nk  =  j->flink! 

t  =  S  - >  flink; 

i  f (t  ==  0)  ( 

s  ->  flink  =  is? 
i  a  =  <*; 
return(t); 
> 
'   t  ->  blink  =  s  ->  blink 
) 

s->flink  =  ia; 
S->blink  =  i  a  -  >  b  1  i  n  k  ; 
i  a  =  s; 
return(t); 
> 


sea  1 e ( ) 
< 

st  rue  t  data  *s  ; 

int  dm,i,j,c»Tiinx,miny; 

int  dx  120)  ,dy  l^ni  ; 


105 


s  = 
d  =  s  i 
fori 


> 

i  f  (m 

i  f  (m 

S  =  db 
for( 


dbuf J  minx  =  miny  =  077777; 


i  =  o; 

minx  = 
miny  = 

s  +  t; 


i  <  np!  i  +  +  )  < 

((c  =  S  ->  x)  <  minx?  c  :  minx); 

((c  =  S  ->  yd)  <  miny?  c  !  miny); 


inx 
m  i  n 
i  ny 
m  i  n 

uf ; 

i  = 
dm 
s  - 
s  - 
s  - 
s  - 
i  f  ( 


> 

s  - 
i  f  ( 

s  - 
s  - 
s  +  * 


>  0) 

x  =  0; 

>  0) 
v  =  o; 

0 ;  i  <  np;  i  ♦  +  )  < 
=  9  ->  yd  K     1  ; 

>  x  =-  minx; 

>  yd  =-  miny; 

>  x  =/  9C 1 x ; 
>x  =♦  biasiJ 
s->x  >=  NRYT  )  < 

print  f ("overflow  pb  #"); 

exi t  ( ); 

>  yd  =/  scly; 
S->yd  >  NSL) 

s->yd  =  NSL-i; 

>  yd  =*  0177  776; 

>  yd  =!  dm; 


i  np ( i  df , buf / nbyte ) 

int  idff*bvjf#nbvtp; 


n  t  I  >  t  f  n  ,  C  » 
t ruct  ( 

char  c 1 »  c2»  c  3. c  u ; 
c  f  f 1  2  0  0 1  ,  *  s  ; 

=  cf; 

if((n  =  read( idf ,c f , nbyte)  )  >  OM 
ford  =  0;  i  <  nbyte/ '4;  i  ♦  +  )  < 

c  =  s  ->  c2  <<  2; 

t  =  c  <  <  io; 

t  =!  s  ->  c  3  <<  6; 

t  = !  s  -  >  c  <j  ; 

s+  ♦ ; 

buf  (i ]  =  t ; 

> 


) 


char  *cbp,  Hoc,  t  Ml  32)  ,cb  110]  ; 


con v ( va 1  ) 

int  va 
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i  nt  a  ; 

i  f  (a  =  val /10) 
con v ( a ) ; 

*cbp++  =  val  2  10  * 


cone (c 1 fc2<n) 
char  «c2; 
i  nt  n,*cL 


I 


l  n  t  i  ,  m  i 

for(i  =  0;  i  <  n;  (++){ 

i  f ( i  ==  0) 

m  =  (*cl  %    07700)  >>  6  ; 
e  1  se 

m  =  *c !  K    0077; 

i  f (m  ==  01?) 

*c2  =  '  0'  ; 

if  (m  >=  01  *,%    m  <=  Oil) 

*e2  =  ' 1 •  +  m  -  l ; 

if  (m  >=  021  &R,  m  <=  051  ) 

*c2  =  'A'  +  m  -  021; 

1  f  f m  >=  OH  I  &R  m  <=  051  ) 

*c2  =  'J'  ♦  m  -0U1; 

i  f  (m  >=  062  &R  m  <=  07  1) 

*c2  =  *S'  +  m  -  062; 

c2  +  * ; 

> 

ret  urn (i  ) ; 


sk i  p (cnt  ) 

I nt  cnt  ; 


< 


i  nt  i  ; 

f  o  r  (  i  =  o;  i  <  cnt;  i  +  I  ) 
write (pdev , pb , 2 ) ; 


clr() 
< 

i  n  t  i  ; 

for(i  =  0;  i  <  132;  i++) 
1 1  f  i )  =  •  •  ; 

1 1  ( l  3 1 1  =  ' n ' ; 
) 


mo v (to» f  rom , n ) 
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char     *  t  o  »  *  f rom  » 
t  nt    n; 

int     \  i 

for(i     =    o;     i     <    n;     i+t) 
* t o  +  +    =     *  f rom** ' 
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